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Abstract 


Since  the  1960's  many  adaptations  to  the  linear  parabolic  model  for  silicon 
oxidation  have  been  proposed.  For  the  purposes  of  process  engineering,  curve  fitting 
procedures  which  are  sometimes  devoid  of  physical  content  but  numerically  precise  are 
employed.  However,  a  truely  unified  physical  model  which  is  in  quantitative  accord  with 
all  known  facts  is  still  lacking.  In  this  review  we  will  discuss  the  "facts"  and  newer 
models . 
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Introduction 


'The  great  technological  importance  of  the  process  by  which  the  surface  of  a  single 
crystal  wafer  of  Si  is  converted  to  a  thin  film  of  amorphous  SiO.  is  derived  from  the 
resultant  electrical  stability  of  the  oxidized  Si  surface.  More  than  twenty  five  years 
ago,  it  was  demonstrated  (1)  that  the  band  gap  states  on  a  Si  surface,  which  arise  as  a 
result  of  the  unsatisfied  bonds,  could  be  reduced  by  about  five  orders  of  magnitude  by 
the  growth  of  adventitious  SiO^ on  the  Si  surface.  Such  a  dramatic  reduction  in  the 
electronically  active  surface  states  via  oxidation  renders  the  Si  surface  useable  for 
high  performance  device  applications.  This  important  finding  has  led  to  the  emergence  of 
the  field  of  planar  integrated  circuits  and  indeed  microelectronics  with  Si  as  the 
preeminent  semiconductor.  Thus,  since  the  early  1960's  there  has  been  a  considerable 
research  and  development  effort,  in  order  to  better  understand  and  improve  the  Si 
oxidation  process. v 


As  a  result  of  the  research  activity,  many  Si  oxidation  models  have  been  proposed. 
Some  of  the  models  have  been  based  on  a  new  fact  or  facts  that  have  been  uncovered  in 
research.  Other  models  have  been  derived  through  analogy  with  similar  materials  that 
have  been  carefully  studied,  and  still  other  models  have  been  proposed  from  curve 
fitting  exercises  where  the  the  concordant  shapes  of  the  model  and  the  oxidation  data  in 
terms  of  SiO^  film  thickness,  L,  versus  oxidation  time,  t,  are  taken  as  confirmatory  of 
the  model.  Most  models  are  combinations  of  these  methods.  It  is  the  purpose  of  this 
paper  to  review  the  Si  oxidation  modelling  effort.  The  next  section  of  this  review 
entitled  "The  Nature  of  The  Oxidation  Process"  will  treat  the  commonly  known  facts  about 
Si  oxidation  such  as  where  the  oxidation  reaction  takes  place  and  characteristics  of  the 
oxide  itself.  Following  this  section  will  be  a  brief  review  of  the  early  models  entitled 
"Early  Models  of  The  Oxidation  Process"  which  includes  a  description  of  the  shape  of  the 
data  and  a  discussion  of  the  typical  experimental  apparatus.  It  will  be  shown  that  many 
aspects  of  the  early  models  are  still  valid.  However,  new  facts  are  constantly  being 
added  to  the  literature  and  these  facts  undoubtedly  lead  to  revised  and  new  models.  The 
next  section  of  this  paper  entitled  "New  Facts  Related  to  the  Oxidation  Process"  will 
review  some  of  the  more  recent  studies  which  have  added  new  facts  for  modelling  and 
following  this  section  will  be  a  section  entitled  "New  Models  For  Si  Oxidation"  which 
will  cover  new  and/or  revised  models.  The  theme  to  be  developed  is  that  while  there  has 
been  considerable  progress  both  in  understanding  the  mechanism  of  Si  oxidation  and  in 
developing  the  process  technology,  there  still  exists  many  important  unanswered 
questions.  At  the  present  time  the  known  facts  are  insufficient  to  unambiguously 
determine  what  is  the  actual  Si  oxidation  mechanism.  However,  in  many  cases  the 
presently  available  information  can  be  used  to  eliminate  some  models  and  this  will  be 
discussed  where  appropriate.  The  final  section  will  summarize  the  present  status  of  Si 
oxidation  modelling  and  suggest  possible  future  directions. 


The  Nature  of  the  Oxidation  Process 


Successful  modelling  of  any  reaction  process  requires  a  complete  knowledge  of  the 
physical  and  chemical  mechanisms  by  which  the  reaction  proceeds.  For  the  specific  case 
of  silicon  oxidation,  an  insulating  film  of  SiO-  is  formed  on  a  Si  surface  in  the 
presence  of  an  oxidizing  ambient.  Since  the  oxide  formed  separates  the  two  reacting 
species,  diffusion  of  one  or  both  of  the  reactants  through  the  oxide  layer  must  occur 
prior  to  the  chemical  reaction  to  form  SiO^.  Therefore,  modelling  the  Si  oxidation 
process  requires  understanding  the  structural  properties  of  the  oxide  in  addition  to  the 
transport  and  chemical  properties  of  the  reactants  since  the  oxide  provides  the  medium 
through  which  the  reaction  process  proceeds.  Figure  1  illustrates  the  Si  oxidation 
process,  showing  the  diffusion  and  reaction  sequence.  It  is  well  established  that  the 
oxidant  is  the  diffusing  species  and  that  the  reaction  to  form  SiO^  occurs  at  the  Si- 
SiC^  interface  (2-10).  However,  the  details  of  the  diffusion  process  and  the  nature  of 
the  oxide  formed  remain  as  valid  research  issues.  Our  current  understanding  of  these  two 
aspects  of  the  oxidation  process  will  be  discussed  below. 

Characterizing  the  diffusion  process  involves  not  only  identifying  which  reactant 
is  diffusing,  but  also  determining  how  the  diffusing  species  interacts  with  the  SiO 
network  and  what  type  of  charge,  if  any,  the  oxidizing  species  possesses.  Early  studies 
indicated  that  oxygen  was  the  diffusing  species  in  dry  oxidation  processes  (3,5)  while 
H^O  was  identified  as  the  diffusing  species  for  steam  ambients  (4).  |gcent  studies 
confirm  and  expand  upon  these  early  conclusions.  Results  from  using  0^  to  trace  the 
oxidation  process  demonstrated  that  the  oxidation  proceeds  predominantly  by  the  long 
range  migration  of  molecular  oxygen  without  exchange  of  atomic  (j>^with  the  oxide  network 
(6,9,10).  The  e^geriments  were  conducted  by  oxidizing  Si  in  an  0^  ambient,  changing 
the  ambient  to  0^  and  determining  which  oxide  regions  contained  isotopic  oxygen.  An 
interesting  result  was  that  although  about  93Z  of  the  isotopic  SiCL  was  located  at  the 
Si-SiO^  interface,  indicating  transport  of  0„  without  exchange  in  the  oxide,  t^gre  was  a 
thin  layer  of  isotopic  SiO_  at  the  outer  surface.  This  outer  layer  containing  0  is 
thought  to  be  related  to  the  motion  of  network  oxygen  atoms  exchanging  near  the  external 
surface  (9,11).  Although  the  effect  of  this  mechanism  on  the  overall  transport  process 
is  small,  its  relative  importance  appears  to  increase  in  the  very  early  stages  of 
oxidation  (9).  Similar  studies  for  oxidation  in  steam  ambients  (7,8)  concluded  that 
transport  in  the  SiC>2  film  is  by  interstitial,  dissolved,  molecular  H20  but  that  a 
reversible  exchange  occurs  between  molecularly  dissolved  H20  and  Si-0  bonds  in  the 
network.  The  charge  state  of  the  diffusing  oxidant  species  has  been  a  source  of 
controversy  for  many  years.  Neutral  0  ,  02  ,  and  02  have  been  suggested  as  possible 
oxidant  species  for  dry  oxygen  ambients.  Reaction  via  a  charged  species  versus  a  neutral 
species  would  significantly  affect  the  oxidation  kinetics  particularly  in  the  early 
stages  of  oxidation  (12).  Evidence  for  oxidation  via  a  charged  species  was  reported  by 
Jorgenson  (4)  in  experiments  which  concluded  that  the  oxidation  process  could  be 
accelerated  or  retarded  by  the  application  of  an  electric  field.  However,  Raleigh  (13) 
disputed  this  argument  claiming  that  the  same  results  would  be  obtained  if  the  electric 
field  had  superimposed  an  electrolytic  reaction  upon  the  normal  oxidation  reaction. 
Recent  results  from  an  experimental  study  of  the  effects  of  electric  fields  on  the 
oxidation  process  by  Modlin  and  Tiller,  demonstrate  that  Jorgensen's  results  are  not 
valid  (14).  Interfering  reactions  from  the  platinum  electrodes  used  in  Jorgensen's  study 
negate  the  conclusions  derived  from  that  investigation.  The  Modlin  and  Tiller  results, 
obtained  by  using  a  corona  discharge  to  create  the  electric  field  and  thereby 
eliminating  metal  contact  with  the  oxide,  concluded  that  the  diffusing  species  is 


neutral  for  oxide  thicknesses  above  300  X.  The  transport  mechanism  for  oxide  growth  up 
to  300  X  has  not  been  clearly  determined.  Thus,  the  best  information  to  date  is  that  for 
dry  oxygen  ambients  the  diffusing  species  is  neutral  0^  which  does  not  interact  with  the 
'xide,  while  for  wet  ambients  the  diffusing  species,  molecular  H^O,  rapidly  exchanges 
with  the  SiO^  network. 

The  physical  and  structural  properties  of  the  oxide  influence  transport  of  the 
oxidizing  species  and  the  chemical  reaction  at  the  interface.  Thus,  the  oxide  structure 
becomes  an  important  factor  in  governing  the  overall  reaction  kinetics.  Understanding 
the  effects  of  oxide  structure  on  the  oxidation  process  is  essential  to  evaluating  the 
kinetic  behavior.  Thermal  oxidation  produces  a  uniform,  continuous  film  of  SiC>2  (1,15) 
on  the  Si  surface.  Electron  diffraction  results  indicate  that  the  SiC>2  film  possesses  an 
amorphous  network  structure  (1),  although  crystallites,  identified  as  cx-crystobal ite , 
have  been  observed  (1,16).  Recent  investigations  related  to  the  oxide  structure  have 
centered  on  several  key  issues:  the  oxide  structure  at  the  Si-SiO  interface,  the 
presence  of  micropores  in  the  oxide,  and  the  effect  of  stress,  introduced  at  low 
temperatures,  on  the  oxide  structure. 

The  bulk  oxide  consists  of  a  three  dimensional  random  network  of  SiO,  tetrahedra 
(17),  but  the  structural  properties  of  the  oxide  located  within  30-40  X  or  the  Si-SiO^ 
interface  are  different  and  have  been  the  subject  of  numerous  investigations  (18-25).  It 
is  well  known  and  not  surprising  that  the  transition  from  crystalline  Si  to  Si02  is  not 
sharp.  Early  studies  (18,19)  indicated  that  a  nonstoichiometric  Si-rich  oxide  existed  at 
the  interface  but  could  not  accurately  determine  the  extent  and  morphology  of  this 
region.  Assimilating  data  from  numerous  TEM  and  spectroscopic  experiments,  the  current 
understanding  of  the  interface  region  is  as  follows  (26).  The  transition  from 
crystalline  Si  to  an  amorphous  oxide  is  rather  abrupt  (about  3  X),  but  there  is  a 
transition  region  of  Si-rich  SiO  which  extends  several  monolayers  into  the  amorphous 
oxide  before  a  stoichiometric  Si02  structure  is  established.  Exact  measurements  of  the 
extent  of  this  region  are  limited  by  interfering  problems  associated  with  the 
measurement  techniques.  Although  the  oxide  attains  a  stoichiometric  Si02  composition 
within  5-10  X  of  the  interface,  there  is  an  additional  transition  region,  extending  to 
approximately  30  X,  which  does  not  possess  the  properties  of  bulk  Si02.  This  region 
appears  to  be  strained,  with  a  distribution  of  tetrahedral  rings  different  from  the 
bulk.  Correlation  of  the  structure  and  extent  of  the  transition  region  with  oxidation 
variables  such  as  temperature  and  pressure  has  not  been  reported.  Determining  such  a 
correlation  may  provide  insight  into  factors  affecting  the  oxidation  reaction  at  the 
interface . 

High  resolution  TEM  studies  have  provided  evidence  for  the  existence  of  voids  or 
"pores"  in  dry  (27,28)  and  wet  (J^)  grown  Si02  films.  These  voids,  ranging  from  5-10  X 
in  diameter  with  a  density  of  10  cm  ,  are  possibly  elongated  in  the  oxide.  Although 
no  direct  evidence  exists  to  indicate  micropores  extending  the  entire  thickness  of  the 
film,  if  they  extend  a  significant  distance  into  the  oxide  they  could  provide  an 
alternate  mode  of  transport  to  the  interface  (29-31).  There  are  also  reports  of 
microchannels  in  the  oxide  which  arise  from  local  ordering  of  the  SiO^  tetrahedra  to 
form  structural  channels  in  the  amorphous  network  (11).  These  microchannels  are  not 
related  to  the  micropore  structures  discussed  above  which  are  gross  structural  defects 
involving  broken  bonds.  No  direct  experimental  evidence  exists  for  microchannels  nor  has 
an  oxidation  model  been  proposed  which  utilizes  these  features. 

The  effect  of  stress  on  the  oxide  structure  has  become  an  issue  of  primary  interest 
as  the  trend  toward  the  use  of  lower  oxidation  temperatures  increases.  There  is  evidence 


that  the  bulk  properties  of  the  oxide  grown  at  temperatures  below  800°C  are  different 
from  those  of  the  high  temperature  oxides  (64).  The  role  of  stress  on  the  oxidation 
process  will  be  discussed  in  detail  in  a  later  section,  but  it  is  important  to  emphasize 
that  the  physical  properties  of  the  oxide  affect  both  the  transport  and  reaction 
processes.  Any  changes  in  the  structural  properties  of  the  oxide  due  to  oxidation 
temperature  or  pressure  will  be  reflected  in  the  oxidation  kinetics. 
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Early  Models  of  the  Oxidation  Process 


As  stated  previously  and  shown  in  Figure  1,  the  oxidation  reaction  proceeds  by  the 
transport  of  oxidant  through  the  oxide  followed  by  reaction  with  silicon  at  the  Si-SiO^ 
interface.  The  overall  rate  of  the  reaction  will  be  governed  by  the  slowest  process. 

Very  early  attempts  to  model  silicon  oxidation  produced  much  debate  over  which 
mathematical  expression  best  represented  the  data.  Numerous  studies  over^a^wide  range  of 
oxidizing  temperatures  and  pressures  used  expressions  of  the  f  rm  L  =  kt  *  (1,2,15,32- 

35),  I,  =  klog(l+at)  (33)  ,  or  L  *  kt  (36)  to  represent  oxide  growth  data  collected  as 
thickness,  L,  versus  time,  t.  These  studies  lacked  a  common  basis  for  comparison  and 
thus  controversy  over  the  description  of  the  oxidation  kinetics  persisted.  However,  most 
studies  in  the  temperature  range  1000°-1300°C  agreed  gh^t  the  data  could  be  reasonably 
well  represented  by  the  parabolic  expression  ,  L  =  kt  *  ,  which  is  consistent  with  a 
diffusion  limited  process.  It  was  recognized,  though,  that  the  data  adhered  to  a 
strictly  parabolic  expression  only  for  oxidation  temperatures  above  1000°C  and  oxide 
thicknesses  greater  than  2000  X  £15).  For  lower  temperatures  and  thin  oxides  a  mixed 
linear-parabolic  expression,  k  L  +  kjL  *  t,  was  required  to  represent  the  data.  This 
mixed  linear-parabolic  form  indicated  that  the  oxidation  was  changing  from  a  purely 
diffusion  limited  reaction  to  one  that  was  also  controlled  by  the  interface  reaction.  In 
a  1965  paper,  which  has  become  a  classic  in  the  field,  Deal  and  Grove  (37)  presented  a 
comprehensive  study  of  oxidation  kinetics  over  a  broad  range  of  temperatures  (700°- 
1200 °C)  and  pressures  (0. 1-1.0  atm)  in  both  wet  and  dry  ambients.  The  model  they 
proposed  to  explain  the  linear-parabolic  nature  of  the  data  has  served  as  the  ba.  ) *»  for 
the  interpretation  of  all  subsequently  collected  Si  oxidation  data.  The  model  in’-*,  es  a 
steady  state  analysis  of  the  physical  and  chemical  processes  known  to  occur  durii. __  the 
reaction.  Figure  2  shows  the  three  processeses  considered  in  the  model:  1)  transport  of 
the  oxidizing  species  across  the  ambient-SiO  interface,  2)  Fickian  diffusion  through 
the  oxide  and  3)  reaction  with  Si  at  the  Si-SiOj  interface.  Under  steady  state 
conditions  these  three  fluxes,  occuring  in  series,  will  be  equal  at  all  times  with 


F .  =  F  *  F 
1.  r2  r3 

The  oxidation  rate  is  given  by  the  equation: 

dL/dt  =  F/a 

where  (1  represents  the  number  of  oxidant  molecules  incorporated  into  a  unit  volume  of 
oxide.  Using  the  flux  equations  presented  in  Figure  2,  the  integrated  form  of  the  rate 
equation  is 


<L2  - 
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where  k  =2DC  /(I  and  k^k^C  /Cl  represent  physico-chemical  parameters  of  the  paction. 
Here  D  ?s  the  effective  diffusion  constant  for  transport  through  the  oxide,  C  is  the 
equilibrium  concentration  of  oxidant  in  the  SiC^,  and  kr  represents  the  reaction  rate 
constant  for  the  conversion  of  Si  to  SiOj.  Transport  across  the  ambi  lt-SiOj  interface 
is  much  faster  than  the  other  two  fluxes  and  thus  the  transport  coefficient,  h,  does  not 
affect  the  reaction  rate  constants  k  or  k^.  The  parabolic  rate  constant,  k  ,  is  related 
to  diffusion  through  the  oxide  while*5 the  linear  rate  constant,  k^ ,  directly^reflects  the 
reaction  at  the  Si  surface.  The  constants  L  and  t  are  offsets  co  the  Deal-Grove  model 
which  account  for  an  initial  region  in  dry  oxidations,  extending  to  about  300  X, 
where  the  model  does  not  accurately  describe  the  experimental  data.  An  offset  is  not 
required  when  wet  ambients  are  used.  Aside  from  the  initial  regime  in  dry  0 ^  ambients, 
the  model  showed  excellent  agreement  for  both  steam  and  dry  oxidation  data.  The 
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oxidation  rate  in  the  initial  regime  is  more  rapid  than  predicted  by  the  Deal-Grove 
model  and  has  been  the  subject  of  numerous  investigations  over  the  past  twenty  years. 

Thus,  although  the  Deal-Grove  model  was  a  significant  contribution  to  the 
interpretation  of  Si  oxidation  data,  it  is  not  complete  and  needs  revisions  particularly 
for  the  initial  stages  of  dry  oxidation  kinetics.  Other  problems  associated  with  the 
model  will  be  discussed  in  the  next  section.  The  importance  of  the  initial  regime 
increases  as  the  oxidation  temperature  is  lowered,  since  at  temperatures  below  800°C  the 
oxide  grown  in  a  reasonable  time  is  entirely  within  the  thickness  of  the  initial  regime. 
The  current  trend  towards  the  use  of  thinner  films  and  lower  processing  temperatures  has 
intensified  the  drive  to  explain  the  oxidation  kinetics  in  this  thickness  regime  (38). 

Before  proceeding  with  the  next  section  it  may  be  useful  to  briefly  describe  the 
procedure  used  for  studying  the  thermal  oxidation  of  silicon.  A  typical  oxidation  is 
performed  in  a  double  walled  quartz  furnace  tube,  as  shown  in  Figure  3,  at  temperatures 
ranging  from  800°-1200°C.  The  furnace  temperature  is  set,  equilibrated  and  calibrated 
prior  to  introducing  the  Si  wafers.  By  flowing  an  inert  gas  such  as  nitrogen  through  the 
outer  jacket  of  the  oxidation  tube  the  diffusion  of  impurities  such  as  Na  and  HO  is 
minimized.  An  inert  ambient,  usually  nitrogen  or  argon,  flows  through  the  inner  cube 
prior  the  oxidation  to  minimize  the  growth  of  any  oxide  during  the  time  it  takes  for  the 
samples  to  reach  the  temperaure  of  the  ambient.  The  cleaning  procedure  for  the  wafers 
before  entry  into  the  oxidation  furnace  is  rigorous  in  order  to  preclude  the  effects  of 
impurities  on  the  oxidation  process.  The  samples  are  then  placed  on  a  quartz  holder, 
entered  into  the  inner  tube  of  the  furnace,  and  pushed  to  the  center  of  the  hot  zone. 
After  allowing  the  samples  to  attain  the  desired  oxidation  temperature  (approximately  15 
minutes),  the  oxidizing  ambient  is  introduced  into  the  inner  tube.  For  dry  oxidations 
this  will  be  pure,  dry  oxygen  while  for  wet  oxidations  the  oxygen  is  first  bubbled 
through  deionized  water.  The  oxidation  is  then  allowed  to  proceed  for  the  desired  time 
interval.  Measurement  of  the  SiOj  film  thickness  is  usually  performed  by  el lipsometry , 
an  optical  technique  based  on  changes  in  the  state  of  polarization  of  an  incident 
monochromatic  light  beam  upon  reflection  from  the  film  covered  Si  surface.  The  advent  of 
the  use  of  in-situ  ellipsometric  studies  where  an  oxidation  furnace  is  mounted  at  the 
focus  of  an  ellipsometer  has  been  an  important  development  in  studying  Si  oxidation 
kinetics  (39-41).  This  in-situ  set  up  enables  the  collection  of  dense  thickness  versus 
time  data  sets  and  therefore  allows  an  accurate  determination  of  the  rate  constants. 


New  Facts  Related  to  the  Oxidation  Process 


The  Deal-Grove  model  (37)  has  repeatedly  provided  an  excellent  fit  to  Si  oxidation 
data  (40-49).  However,  in  addition  to  the  unexplained  initial  regime,  detailed  studies  of 
various  aspects  of  the  oxidation  process  have  uncovered  new  facts  which  are  either 
inconsistent  with  the  Deal-Grove  model  or  not  accounted  for  in  the  model.  The  use  of  in- 
situ  el  1 ipsometric  studies  has  been  instrumental  in  revealing  details  of  the  oxidation 
mechanism.  These  studies  and  others  support  the  Deal-Grove  model  as  a  basis  for  analyzing 
Si  oxidation  data,  but  indicate  that  the  oxidation  process  is  more  complex  than 
originally  proposed.  In  this  section  the  properties  of  the  rate  constants  and  overall 
reaction  rate  predicted  from  the  model  will  be  discussed  and  compared  with  new  facts 
derived  from  observed  experimental  behavior.  Attempts  to  incorporate  the  new  facts  into 
models  of  the  oxidation  process  will  be  discussed  in  the  next  section. 

In  the  Deal-Grove  model  the  parabolic  rate  constant,  k  =  2DC*/(1,  reflects  the 
diffusional  transport  of  oxidant  through  the  oxide.  If  k  corresponds  to  an  elementary, 
thermally  activated  process  obeying  Boltzman  statistics .^its  temperature  dependence 
should  be  described  by  an  Arrhenius  equation  of  the  form 

k  =  kQexp(-E  /kgT) 

where  kfi  is  the  Boltzman  constant,  T  *is  tlFe  oxidation  temperature,  and  E  is  the  single 
activation  energy  for  the  process  which  in  this  case  is  diffusion  through  the  oxide.  For 
a  kinetical ly  simple  rate  process  described  by  a  single  rate  constant,  say  k  ,  the 
Arrhenius  plot,  ln(k  )  versus  1/T,  will  be  a  straight  line  with  a  slope  equa?  ro  the 
activation  energy.  Dial  and  Grove  showed  that  the  k  values  derived  from  their 
experimental  results  exhibited  Arrhenius  behavior  wfth  activation  energies  determined  for 
wet  and  dry  oxidizing  ambients  that  were  similar  to  those  measured  for  the  diffusion  of 
water  and  oxygen,  respectively,  in  fused  silica.  They  also  reported  a  linear  pressure 
dependence  for  k  over  the  range  0. 1-1.0  atm  from  an  analysis  of  literature  data.  This 
linear  pressure  iependence  is  concordant  with  the  model  since  k  is  directly  proportional 
to  C  ,  the  equilibrium  concentration  of  oxidant  in  the  oxide.  According  to  Henry's  law, 
in  the  absence  of  dissociation,  the  equilibrium  concentration,  and  thus  k  ,  will  be 
linearly  related  to  pressure.  The  pressure  dependence  reported  for  k  is  therefore 
consistent  with  diffusion  via  a  molecular  species.  Deal  and  Grove  concluded  that  the 
model  successfully  predicted  the  temperature  and  pressure  behavior  of  k  .  Further  studies 
have  confirmed  the  linear  pressure  dependence  of  k  at  temperatures  ranging  from  870°- 
1200°C  (41,50).  There  is  some  question  about  the  pFessure  dependence  at  lower 
temperatures,  but  the  limited  growth  rate  makes  it  difficult  to  determine  rate  constants 
accurately.  The  behavior  of  k  with  temperature  has  been  carefully  studied  and  it  is  now 
apparent  that  k  does  not  exhFbit  the  simple  Arrhenius  behavior  originally  reported  in 
the  Deal-Grove  itudy  (48,49,51).  While  the  activation  energies  reported  in  the  Deal-Grove 
study  are  valid  for  oxidation  temperatures  above  J000°C,  these  new  low  temperature 
studies  show  a  change  in  the  slope  of  the  Arrhenius  plot  occuring  in  the  range  900°- 
1000°C  (48,52).  At  lower  temperatures  the  slope  of  the  curve  (i.e.  the  activation  energy) 
increases  and  thus  the  shape  of  the  plot  is  slightly  concave  downwards.  This  curvature  is 
consistent  with  the  shape  of  an  Arrhenius  plot  for  a  complex  rate  constant  describing  two 
processes  with  different  activation  energies  occuring  in  parallel.  However,  it  does  not 
prove  that  two  transport  mechanisms  are  involved,  k  was  observed  to  be  orientation 
independent  at  temperatures  above  1000°C  (42)  as  anFicipated  for  diffusion  through  an 
amorphous  oxide.  However,  at  lower  temperatures  k  varies  with  Si  orientation  for  both 
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wet  and  dry  oxidizing  ambients  (47,49).  Thus  it  appears  that  the  diffusion  process  or  the 
structure  of  the  oxide  is  influenced  by  the  surface  structure  of  the  Si. 

The  linear  rate  constant,  k^  =  k  C  /(I,  should  be  a  direct  reflection  of  the 
reaction  at  the  surface.  Any  changes  in  the  reaction  rate  constant,  k  ,  will  directly 
affect  k^.  Deal  and  Grove  reported  that  k.,  in  addition  to  k  ,  exhibited  Arrhenius 
behavior.  The  activation  energy  determined  for  k^  was  similaF  to  the  Si-Si  bond  energy 
indicating  that  this  bond  breaking  step  was  rate  limiting.  The  Deal-Grove  igodel  also 
predicts  a  linear  pressure  dependence  for  k^  since  it  is  proportional  to  C  .  New  findings 
present  a  more  complex  picture  of  the  linear  rate  constant:  the  Arrhenius  plot  for  k  is 
not  linear  (45,48,49)  and  k.  does  not  exhibit  a  linear  pressure  dependence  (41,46,49). 

The  Arrhenius  plot  for  k^,  tike  the  plot  for  k  ,  exhibits  a  break  in  the  slope  occuring 
between  900°C  and  1000°C,  but  the  direction  of ^curvature  is  opposite  to  that  for  k  .  A 
higher  activation  energy  is  observed  for  k^  at  temperatures  above  1000°C  than  for  ¥ower 
temperatures  (45,48,49).  This  shape,  concave  upward,  parallels  the  shape  of  an  Arrhenius 
plot  for  a  complex  rate  constant  representing  processes  with  different  activation 
energies  occuring  in  series  (52).  Again,  the  shape  of  the  curve  is  not  sufficient 
evidence  to  prove  this  type  of  complex  mechanism  exists.  The  pressure  dependence  of  k^ , 

P  ,  varies  with  temperature  from  n=0.59  at  low  temperatures  to  n=0.83  at  high 
temperatures  (45).  The  exponent,  n,  varies  with  orientation  for  a  fixed  temperature, 
which  also  indicates  that  k^  does  not  represent  a  simple  reaction  step  (45).  This  complex 
pressure  dependence  has  been  explained  in  terms  of  reaction  by  both  molecular  and  atomic 
species  at  the  interface.  Additionally,  k^  exhibits  a  distinct  orientation  dependence  as 
a  function  of  temperature  (45,47,49)  at  1  atm.  It  has  been  proposed  (53)  that  the 
relative  ordering  of  k  values  is  related  to  the  number  of  Si-Si  bonds  on  the  surface  of 
a  particular  plane  with  (110)  >  (111)  >  (100).  However,  this  relative  ordering  of  k 
values  only  holds  for  temperatures  above  900°C  (47,49)  and  thus  other  factors  are  also 
involved . 


Interestingly,  both  k  and  k  exhibit  breaks  in  their  Arrhenius  plots  in  the 
temperature  range  900°-1008°C,  although  the  direction  of  curvature  is  opposite  for  the 
two  cases.  This  change  in  slope  reflects  a  change  in  the  oxidation  mechanism  at  low 
temperatures.  As  previously  mentioned,  it  could  indicate  the  presence  of  additional 
processes  occuring  in  parallel  or  in  series,  respectively,  with  the  diffusion  and 
reaction  fluxes  stated  in  the  Deal-Grove  model  (52).  If  additional  fluxes  are  not 
present,  then  the  curvature  could  reflect  a  thermally  controlled  property  influencing  the 
transport  and  reaction  steps  described  in  the  model.  Evidence  for  a  temperature  dependent 
stress  affecting  the  kinetic  behavior  is  available  and  will  now  be  considered.  Upon 
cooling  from  the  oxidation  temperature,  a  stress  develops  in  the  oxide  due  to  differences 
in  the  thermal  coefficients  of  expansion  for  Si  and  SiO_  (55).  The  temperature  reported 
for  relaxation  of  this  thermal  stress  through  viscous  flow  is  around  970°C  (55).  This 
temperature  lies  precisely  within  the  range  where  a  break  occurs  in  the  Arrhenius  plots 
for  both  kj  and  k  ,  and  thus  indicates  viscous  flow  may  play  an  important  role  in  the 
oxidation  process?  Further  studies  revealed  th^t  in  addition  to  the  thermal  stress  an 
intrinsic  stress,  on  the  order  of  10~  dynes/cm  ,  exists  in  oxides  grown  in  both  wet  (56) 
and  dry  (56,57)  ambients  at  temperatures  below  950°C.  The  intrinsic  stress  is  a  stress 
component  which  exists  at  the  oxidizing  temperature  where  the  thermal  stress  component 
zero.  This  intrinsic  stress  probably  arises  because  the  oxide  lacks  the  ability  for 
viscous  flow  to  occur  at  these  temperatures  (57),  and  is  not  observed  at  oxidizing 
temperatures  above  975°C  (56).  Thus,  a  new  parameter,  intrinsic  stress,  is  introduced 
into  the  oxidation  process  in  the  temperature  range  950°-975°C  and  is  present  at  all 
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lower  oxidation  temperatures.  Although  no  direct  correlation  has  been  made  between  stress 
and  the  rate  constants,  the  effects  of  stress  on  the  oxidation  reaction  have  been 
considered  (31,58-61).  It  has  also  been  reported  that  the  refractive  index  of  thermal 
SiC^  increases  as  the  oxidation  temperature  is  lowered  (57,62,63).  This  indicates  a 
structural  change  in  the  oxide  as  the  temperature  is  lowered.  Annealing  the  oxide  in  an 
inert  ambient  at  temperatures  where  viscous  flow  can  occur  allows  the  oxide  to  relax  to  a 
refractive  index  characterist ic  of  high  temperature  thermal  oxides  (57).  The  higher 
refractive  index  was  originally  attributed  to  an  increased  oxide  density  (57,62,63),  but 
very  recent  results  indicate  the  refractive  index  change  may  be  caused  by  an  incomplete 
oxidation  leaving  a  Si-rich  oxide  (64).  Thus,  as  discussed  previously,  the  nature  of  the 
oxidation  process  depends  on  the  properties  of  the  oxide.  The  presence  of  stress  and 
altered  structural  properties  introduced  at  low  temperatures  could  influence  transport 
and  the  reaction.  Their  effect  on  the  kinetic  mechanism  must  be  determined  in  order  to 
develop  a  comprehensive  model  for  the  oxidation  process. 

Other  observations  have  been  made  about  the  overall  reaction  kinetics  in  terms  of 
the  Deal-Grove  model.  Determining  the  pressure  dependence  of  k  and  k^  individually  at 
low  pressures  is  difficult  since  the  limited  oxide  growth  precfudes  extracting  accurate 
rate  constants.  Therefore,  some  studies  have  investigated  the  pressure  dependence  of  the 
overall  reaction  rate  (41,46,49).  According  to  the  Deal-Grove  model  the  oxidation  rate 
should  be  linearly  proportional  to  pressure  over  all  thickness  ranges.  However,  while  the 
pressure  dependence  of  dL/dt  is  found  to  be  linear  for  oxides  greater  than  300  X  (41), 
for  ox^d^s  lesg  ghan  300  X  various  studies  have  reported  a  pressure  dependence  ranging 
from  P  to  P  (41,46,49).  Another  interesting  observation  of  Si  oxidation  involves 
the  relative  rates  at  which  the  various  crystal  orientations  oxidize  as  a  function  of 
oxidation  temperature,  pressure  and  time.  For  all  temperatures  in  the  range  800°-1000°C 
the  (100)  Si  orientation  has  the  slowest  oxidation  rate  (47,65)  at  1  atm.  However,  there 
is  a  crossover  in  the  relative  rate  of  oxidation  between  the  (111)  and  (110)  surfaces  at 
each  temperature  in  this  range  (65).  In  all  cases  the  (110)  orientation  has  the  greatest 
initial  oxidation  rate,  but  at  a  thickness  characteristic  of  the  oxidation  temperature 
the  oxidation  rate  on  the  (111)  surface  becomes  greater  than  on  the  (110)  surface  (65). 
The  relative  rates  at  very  low  pressures  as  a  function  of  temperature  is  quite  complex 
(65,66).  Further  research  on  the  relative  oxidation  rates  involving  additional 
orientations  is  presently  underway  in  our  laboratory. 

Information  on  the  effects  of  impurities  on  the  oxidation  process  has  been  collected 
and  may  provide  some  insight  into  the  oxidation  mechanism.  During  the  evolution  of 
integrated  circuit,  IC,  technology  the  emphasis  on  oxidation  impurity  effects  has 
changed.  The  earliest  studies  were  aimed  at  devising  processing  methods  that  would 
eliminate  those  impurities  that  degrade  IC  device  performance.  The  most  notable  impurity 
was  sodium.  Na  can  diffuse  quite  rapidly  in  SiO^  and  with  its  prodigous  natural 
abundance  is  a  major  Si-SiC^  interfacial  impurity  giving  rise  to  large  amounts  of 
positive  charge  at  the  Si  surface.  In  terms  of  oxidation  kinetics,  Na  has  been  found  to 
enhance  the  Si  oxidation  rate(67).  The  details  of  the  kinetics  are  not  well  known. 

H^O,  as  both  an  impurity  and  as  an  oxidant  has  received  considerable 
attention(47,69).  As  an  oxidant,  HjO  is  a  more  virulent  oxidant  for  Si  than  is  0_.  Hence 
^0  is  used  for  processing  when  thick  oxides  are  required  so  as  to  reduce  both  tne 
process  time  and  temperature.  However,  the  use  of  HjO  gives  rise  to  a  large  number  of 
bulk  oxide  charge  trapping  centers  which,  when  charged,  also  alter  the  Si  surface 
potential  as  does  Na  at  the  interface  and  affects  the  electric  reliability  of  the 
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device(68).  H-0  oxidation  kinetics  have  been  studied  in  greater  depth  than  for  Na+.  It 
was  found  that  even  traces  of  H^O  in  C>2  increase  the  rate  of  oxidation(47 ,69) .  Both  k} 
and  k  are  enhanced.  For  the  overall  rate  of  oxidation  the  amount  of  enhancement  was 
foundpto  be  greater  than  for  0^  and  H_0  present  as  independent  oxidants  ie.,  there  is  a 
cooperative  effect(69).  Essentially,  this  arises  from  the  fact  that  1^0  reacts  with  the 
Si(>2  network  forming  SiOH  species  and  thereby  terminating  the  Si-O-Si  chains .  This 
loosened  network  enables  rapid  diffusion  of  oxidant  species  whether  the  species  are 
related  to  either  C>2  or  l^O.  has  also  been  found  to  seriously  modify  the  transport  of 

C>2  in  other  ways.  Several  complimentary  studies  (6-8,70)  demonstrate  that  pure  0 
diffuses  through  the  interstices  in  the  Si02  network.  For  the  purpose  of  explaining  Si 
oxidation  there  is  no  exchange  of  0  with  the  SiO„  network.  However,  with  even  traces  of 
H2O  present  0  exchanges  rapidly  with  the  Si02  network,  hence  dramatically  altering  the 
transport  of  oxidant  and  the  oxidation  mechanism. 

Other  impurities  added  to  the  oxidizing  ambient  have  been  considered  such  as  HCl, 

CI2  and  organic  chlorine  containing  compounds(71-73).  The  literature  attests  to  the  fact 
that  oxide  charge  levels  are  reduced  and  oxide  reliability  is  enhanced  through  these 
additions.  The  effects  on  Si  oxidation  kinetics  are  about  the  same,  ie.  all  of  these 
impurities  increase  the  rate  of  oxidation.  The  explanation  is  likely  the  same  in  all  of 
these  cases.  The  impurities  all  react  with  the  Si02  network  thus  enabling  more  rapid 
transport  of  oxidant  to  the  Si-SiC>2  interface  where  oxidation  takes  place. 

Impurities  or  dopants  added  to  the  Si  substrate  also  alter  the  oxidation  kinetics. 
The  most  common  Si  dopants  are  n-type  dopants  such  as  P  or  As  and  p-type  dopants  such  as 
B.  Both  kinds  of  dopants  have  been  found  to  alter  the  Si  oxidation  kinetics  but  with 
different  temperature  dependencies (48 , 74 ) .  At  oxidation  temperatures  of  less  than  1000°C, 
the  n-type  dopants  enhance  the  overall  oxidation  rate  while  at  higher  temperatures  B 
doping  results  in  a  higher  rate.  It  appears  that  both  k^  and  k  are  affected.  The  k 
effects  are  likely  due  to  alterations  of  the  SiC^  network  as  aPresult  of  bonding  dePects. 
The  dopants  B,  P,  and  As  are  network  formers  but  the  resultant  networks  are  different  in 
structural  arrangement  and/or  number  of  bonds  to  0(48).  B  forms  a  planar  trigonal 
arrangement  of  0's  while  P  and  As  form  four  bonds  to  0  as  does  Si  but  the  extra  pair  of 
non-bonding  electrons  distorts  the  tetrahedron.  Such  defects  are  likely  to  enhance 
diffusion  of  the  oxidant.  The  interfacial  oxidation  kinetics  are  even  more  complex.  One 
model (75)  utilizes  the  production  of  point  defects  at  the  interface  to  explain  the 
enhanced  rate  for  P  doping.  The  Si  band  gap  modif icat ions  were  shown  to  increase  the 
vacancy  production  and  this  scaled  with  the  oxidation  rate  enhancement.  Another  new  model 
for  surface  rate  enhancement  has  been  proposed ( 76 ) .  In  this  model  the  availabi  lity  of 
electronic  carriers  is  the  important  property.  For  the  case  of  sub-monolayer  gold 
coverage  on  Si,  it  was  shown  that  the  unannealed  surface  was  metallic  while  annealing 
ordered  the  surface  and  the  ordering  gave  rise  to  a  semiconductor-like  surface.  The 
metallic  surface  oxidized  considerably  more  rapidly  than  the  semiconductor  surface 
presumably  due  to  the  availability  of  electrons  for  the  formation  of  bonds. 

Interestingly,  the  oxidation  of  metal  silicides  also  shows  similar  behavior ( 77 ) .  When  a 
metal  silicide  film  is  grown  directly  on  Si,  the  siiicide  will  oxidize  to  form  SiC>2  only 
with  the  preservation  of  the  silicide  film.  In  effect  there  is  a  net  flow  of  Si  to  the 
si 1 icide-oxidant  interface.  It  was  discovered  that  the  metallic-like  silicides  oxidized 
much  more  rapidly  than  the  semiconductor-like  silicides  and  these  latter  silicides 
oxidized  with  similar  rates  as  Si  itself.  It  seems  that  several  experimental  studies  have 
shown  the  importance  of  the  electrical  nature  of  the  oxidizing  surface.  Such  effects  have 
largely  been  ignored  in  many  oxidation  models. 


It  is  apparent  that  a  considerable  amount  of  information  has  been  collected  in  the 
past  20  years  giving  more  insight  into  the  details  of  the  oxidation  process.  These 
investigations  have  not  only  determined  complex  behavior  in  the  rate  constants  but  have 
also  recognized  the  influence  of  additional  physical  phenomena  introduced  at  low 
temperatures  on  the  oxidation  reaction.  Any  model  which  provides  a  complete  description 
of  the  oxidation  process  must  be  able  to  account  for  the  observed  experimental  facts.  It 
is  clear  that  new  models  or  revisions  to  the  Deal-Grove  model  are  necessary  to  provide  a 
comprehensive  description  of  the  oxidation  mechanism. 


New  Models  For  Si  Oxidation 


Structural  Model 


Considering  the  recent  emphasis  on  the  effects  of  intrinsic  stress  on  Si  oxidation, 
Tiller  (78)  proposed  a  model  that  minimizes  the  stress  and  hence  the  free  energy  required 
for  the  oxidation  process.  Essentially  Tiller  considers  oxidation  to  be  a  two  step 
process.  In  the  first  step  Si  is  converted  to  ac-cristobalite  plus  Si  interstitials  in 
both  Si  and  SiO^.  The  ac-cristobalite  lattice  fits  best  on  a  Si  lattice  except  for  the 
presence  of  a  few  Si  atoms,  and  thus  yields  the  smallest  molar  volume  change  during 
oxidation.  This  first  step  minimizes  the  strain  and  hence  the  stress.  In  the  second  step 
the  Si  interstitials  in  the  SiO^  are  oxidized.  The  volume  expansion  inside  the 
crystalline  ac-cristobalite  lattice  causes  distortion  of  the  lattice  and  eventual 
amorphization.  One  implication  of  this  model  is  that  a  layer  of  crystalline  SiO^  should 
always  be  found  at  the  Si-SiO ^  interface.  While  there  is  no  overwhelming  evidence  for 
this,  Aguis  et  al(79)  have  reported  that  under  low  oxidant  pressures  a  crystalline 
modification  of  SiO  is  observed.  Additionally,  virtually  all  the  modern  surface 
spectroscopies (see  ror  example  refs.  18,19,26)  agree  that  the  interfacial  oxide  is  a 
suboxide  of  SiO^.  Some  recent  work  in  our  laboratory  (64)  on  the  structure  of  low 
temperature  grown  thermal  oxide  on  Si  shows  that  suboxides  are  formed  in  greater  amounts 
when  the  stress  is  highest,  ie.  at  the  lowest  oxidation  temperatures.  We  propose  that  the 
oxidation  system  accomodates  to  the  accumulation  of  intrinsic  stress  which  arises  from 
the  large  molar  volume  change  by  rejecting  complete  oxidation  and  thereby  reducing 
somewhat  the  large  molar  volume  change.  Mott(80-82)  has  considered  that  if  oxidation  of 
Si  proceeds  at  kink  sites  in  the  Si  surface,  then  the  large  volume  requirement  is 
diminished  and  allows  oxidation  to  proceed.  The  applicabil ity  of  this  idea  has  not  been 
determined,  but  it  is  worthwhile  to  consider  since  this  mechanism  may  result  in  small 
strains  that  diminish  with  increasing  temperature.  This  model  may  therefore  provide  an 
alternate  explanation  to  the  viscous  flow  model  which  will  be  introduced  in  the  section 
on  stress  models. 


Interface  Kinetics  Models 


Twd  interface  kinetics  models  have  received  at  tent ion( 34 , 83 ) .  The  earlier  Ghez  and 
van  der  Meulen  model (54)  is  based  on  the  reaction  between  Si  and  oxidant  at  the  Si-SiC^ 
interface  involving  both  0.  and  0.  The  transported  oxidant  species  is  neutral  0?.  This 
model  yields  the  experimentally  determined  pressure  dependence  for  k^(45)  and  shows  that 
oxidation  at  the  higher  temperatures  involves  more  0^  while  at  lower  temperatures  0  is 
the  dominant  oxidant.  At  the  interface  there  is  an  equilibration  reaction  to  produce  0 
from  the  arriving  0^.  A  subset  of  this  model  was  proposed  by  Blanc(83)  where  the  initial 
oxidation  regime,  ie .  the  L  ,  t  regime,  could  be  analytically  accounted  for  by 
considering  the  reaction  at  the  Si-SiO^  interface  to  involve  only  0.  However  successful 
this  model  is  at  explaining  the  initial  regime,  it  contains  neither  the  correct  pressure 
dependence  for  k.  nor  the  correct  temperature  dependence.  On  the  temperature  dependence, 
it  was  found  that  the  Arrhenius  plot  for  the  k.  displayed  a  concave  upwards  shape  which 
is  indicative  of  a  parallel  rate  process(52).  This  is  in  agreement  with  the  Ghez,  van  der 


Models  Based  On  Charged  Species 


Several  models  have  been  proposed  which  incorporate  the  effects  of  charged  species 
on  the  oxidation  kinetics.  These  models  are  all  based  on  the  assumption  that  the 
oxidizing  species  possesses  a  negative  charge.  This  assumption  arises  from  results 
reported  by  Jorgensen  (4)  claiming  that  the  oxidation  process  could  be  enhanced  or 
retarded  by  varying  the  direction  of  an  electric  field  applied  to  the  oxide.  Recent 
results  by  Modlin  and  Tiller  (14),  discussed  in  the  section  on  the  nature  of  the 
oxidation  process,  demonstrate  that  Jorgensen's  results  are  not  valid,  and  concluded  that 
the  diffusing  species  is  neutral  for  oxide  thicknesses  greater  than  300  X.  Although  these 
results  seem  to  dispute  any  models  based  on  oxidation  by  a  charged  species,  they  are  not 
conclusive  for  oxide  thicknesses  below  300  X.  Thus  it  is  useful  to  present  a  brief  review 
of  these  models  since  they  represent  some  thoughts  on  the  source  of  the  rapid  initial 
regime . 

The  initial  rapid  oxidation  rate  observed  for  dry  ambients  was  recognized  by  Deal 
and  Grove  when  they  presented  their  model  (37).  To  account  for  this  behavior  they  applied 
a  theory  developed  by  Cabrera  and  Mott  (12)  to  explain  the  initial  rapid  oxidation  rate 
observed  for  metals.  The  theory  predicts  that  space  charge  effects  will  control  the 
oxidation  until  a  thickness  is  reached  where  electron  transport  through  the  oxide  due  to 
tunneling  or  thermionic  emission  is  no  longer  possible.  This  thickness,  given  by  the 
Debye  length,  is  approximately  150  X  for  1000°C  dry  oxidations.  Deal  and  Grove  thus 
concluded  that  the  Mott-Cabrera  theory  provided  a  possible  explanation  for  the 
anomolously  rapid  rate.  However,  it  must  be  remembered  that  this  theory  was  developed  for 
the  oxidation  of  metals  where  a  readily  available  source  of  electrons  exists  at  the  metal 
surface  and  for  a  conducting  oxide  which  could  transport  charge  and  maintain  overall 
charge  neutrality.  A  later  model  proposed  by  Grove  (84)  attributed  the  enhanced  rate  to 
field  effects  created  by  the  formation  of  an  0-  -  hole  pair  when  the  oxygen  molecule 
enters  the  oxide.  The  greater  mobility  of  the  Role  creates  an  additional  "pull"  on  the 
ion  and  thus  enhances  the  diffusion.  Other  models  have  been  proposed  which  assume  the 
charged  oxidizing  species  is  influenced  by  the  presence  of  oxidation  induced  charges  in 
the  oxide  (86-89).  Positive  interfacial  electronic  charge,  so-called  fixed  oxide  charge, 
Q^,  is  known  to  exist  at  the  Si-SiO^  interface  (85).  One  model  (87)  proposes  that  this 
fixed  oxide  charge  establishes  an  electric  field  within  150-200  X  of  the  interface 
region,  and  thus  enhances  the  diffusion  of  the  negatively  charged  oxidizing  species. 
Incorporating  an  electric  field  assisted  diffusion  process  into  the  Deal-Grove  model,  a 
linear-parabolic  expression  is  derived  with  rate  constants  that  are  exponent ial ly 
dependent  on  the  charge  distribution.  Two  fitting  parameters  related  to  the  charge 
distribution  enable  a  fit  to  the  data.  Models  by  Hamasaki  (88)  and  Naito  et.  al.  (89) 
argue  that  the  electric  field  established  at  the  interface  will  inhibit  the  negatively 
charged  oxidant  species  from  reaching  the  interface.  The  most  recent  model  (86)  considers 
that  the  initially  large  number  of  Si  bonds  at  the  bare  Si  surface  is  reduced  by  the 
presence  of  via  the  reduction  in  the  number  of  holes  by  band  bending.  Neither  the 
origin  of  Qf  nor  the  reasons  why  broken  bonds  are  crucial  are  treated  in  the  model.  From 
an  assumed  distribution  of  charge  and  the  1-D  solution  to  Poisson's  equation,  a  model  is 
derived  which  when  fitted  to  the  oxidation  data  yields  physically  acceptable  values  for 
two  parameters,  the  fixed  oxide  charge  and  the  extent  of  the  initially  fast  oxidation 


regime.  None  of  these  models  has  been  tested  either  theoretically  or  experimentally 
against  the  myriad  of  information  about  such  as  its  orientation  and  ambient  dependence 
and  the  commensurate  effects  on  oxidation.  Therefore,  at  the  present  such  models  should 
be  treated  as  interesting  curve  fitting  exercises. 


Stress  Models 


Since  the  discovery  of  a  considerable  intrinsic  compressive  stress  in  oxides  grown 
at  temperatures  below  950°C  (56)  numerous  models  have  been  proposed  to  explain  the  effect 
of  stress  on  the  oxidation  reaction.  The  controversy  over  the  role  of  stress  centers  on 
whether  stress  affects  the  diffusion  process  or  the  reaction  at  the  interface.  Models 
have  been  proposed  for  both  cases,  but  before  discussing  these  models  it  is  useful  to 
present  a  model  which  describes  the  origin  of  the  intrinsic  stress  (57).  The  conversion 
of  Si  to  SiO.  involves  a  2.2  fold  increase  in  molar  volume.  Referring  to  Figure  4  it  is 
apparent  that  the  oxide  is  constrained  to  the  surface  of  the  wafer,  and  thus  the 
additional  volume  required  for  the  formation  of  SiO^  must  be  obtained  in  the  direction 
normal  to  the  surface.  At  high  temperatures,  where  viscous  flow  occurs,  this  additional 
volume  is  readily  attained  since  the  oxide  can  easily  flow  in  this  direction.  However,  at 
temperatures  below  950°C,  where  viscous  flow  is  not  observed  (55),  the  oxide  cannot 
attain  the  required  molar  volume,  and  thus  an  intrinsic  stress  develops.  Whether  this 
intrinsic  stress  predominantly  affects  the  kinetic  behavior  of  the  diffusion  process  or 
the  reaction  to  form  SiO^  has  been  a  topic  of  much  debate.  One  model  described  below 
emphasizes  the  effect  of  stress  on  the  reaction  at  the  interface,  while  two  other  models 
consider  the  effect  of  stress  on  diffusion  but  reach  different  conclusions. 


A  model  proposed  by  Irene  (31)  introduces  the  effect  of  stress  through  the  linear 
rate  constant.  This  model  contends  that  the  lack  of  viscous  flow  at  low  oxidation 
temperatures  inhibits  the  reaction  at  the  interface  by  reducing  the  number  of  Si  atoms 
that  are  available  for  reaction.  The  formation  of  SiO.  without  the  ability  to  attain  the 
required  molar  volume  in  the  normal  (z)  direction  will  screen  other  Si  atoms  on  the 
surface  from  reacting.  A  modified  expression  f<y  Jfhe  reaction  flux  is  given  by 
*  F3*  krC  Cg. 

where  Cg^  represents  the  effective  concentration  or  Si  atoms  available  for  reaction.  This 
concentration  is  then  controlled  by  the  rate  at  which  SiO.  can  viscously  relax  and 
thereby  reduce  the  screening  effect.  Treating  the  oxide  as  a  Maxwell  solid,  the 


expression  derived  for  C 


Cc.  -  k'Ce.o-  !\) 
■face,  tt  ^.s  fflfe  vi 


iscosity  of  the  SiO„,  and  C, 


where  o-  is  the  stress  at  the  interface,  is  fne  viscosity  of  the  SiO.,  and  C  .  i 
the  actual  concentration  of  Si  atoms  on  the  surface.  Incorporating  this  revised  flux 
expression  into  the  Deal-Grove  model,  the  integrated  rate  equation  retains  a  linear- 
parabolic  form  but  the  linear  rate  constant  tlbecoige8 


k  C„.C  cr 
Si  xv 


The  revised  linear  rate  constant  introduces  an  orientation  dependence  through  C  .  and 
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o-^  ,  and  predicts  a  complex  temperature  behavior  through  o-^  and  ft.  The  behavior 
predicted  by  this  model  has  not  been  confirmed  because  accurate  values  for  o-  and  ft  as 
a  function  of  orientation  and  temperature  have  not  been  reported.  However,  the^  model  is 
based  on  physical  principles,  and  thus  may  provide  an  qualitative  if  not  a  quantitative 
description  of  the  effects  of  stress  on  the  interface  reaction.  This  model  has  been  used 
to  explain  the  crossover  in  oxidation  rate  observed  between  the  (110)  and  (111) 
orientations  at  temperatures  between  800°C  and  1000°C.  It  was  suggested  (65)  that  the 
initial  oxidation  rate  is  governed  by  a  surface  related  property  such  as  the  number  of  Si- 
Si  bonds.  Beyond  the  crossover  point,  the  oxide  is  thick  enough  to  exert  a  force  on  the 
interface,  arising  from  the  intrinsic  stress,  which  is  sufficient  to  alter  the  reaction 
kinetics.  Thus,  the  reaction  beyond  the  crossover  is  influenced  by  the  stress  state  of 
the  oxide. 

Two  models  have  been  proposed  which  claim  that  stress  affects  the  diffusion  process 
rather  than  the  interface  reaction,  but  the  functional  forms  predicted  by  the  two  models 
for  L  versus  t  are  quite  different.  The  model  proposed  by  Doremus  (58)  asserts  that  the 
linear  part  of  the  linear-parabolic  expression  does  not  represent  the  reaction  at  the 
interface  but  instead  arises  from  strain  in  the  oxide.  The  model  is  based  on  the 
assumption  that  the  only  parameter  controlling  the  oxidation  rate  is  the  diffusion  of 
oxidant  to  the  interface.  The  model  assumes  an  exponentially  decaying  diffusion  constant 
as  a  function  of  strain  in  the  oxide: 

D  -  D  exp(-k6o(l  -  J)) 

where  6  represents  the  maximum  strain  located  at  the  interface ,  DL  is  the  diffusion 
cons  tan?  in  the  outer  unstrained  oxide,  L  is  the  total  oxide  thickness  and  k  is  a  fitting 
parameter.  Using  this  expression  for  the  diffusion  constant,  a  linear-parabolic  equation 
is  derived  for  thickness  as  a  function  of  time.  However,  the  linear  term  is  introduced 
only  by  assuming  that  €  ,  the  maximum  strain  located  at  the  Si-SiO?  interface,  is 
inversely  proportional  ?o  the  oxide  thickness.  Thus,  this  model  only  holds  if  this  is  the 
correct  form  for  strain  at  the  interface  as  a  function  of  thickness.  This  model  does  not 
account  for  the  observed  pressure  dependence  of  the  linear  rate  constant  and  more 
importantly  fails  to  account  for  the  fact  that  a  linear  term  is  required  to  fit  the 
experimental  data  even  at  temperatures  above  950 °C  where  no  strain  in  the  oxide  is 
observed.  The  other  model  which  considers  the  effect  of  stress  on  the  diffusion  process 
introduces  stress  directly  into  the  parabolic  rate  constant  of  the  Deal-Grove  model.  This 
model  (60),  proposed  by  Fargeix  and  Ghibaudo,  is  based  on  a  curve  fitting  procedure 
applied  to  plots  of  dt/dL,  the  inverse  growth  rate,  versus  L.  A  typical  plot  of  dt/dL 
versus  L  for  experimental  data  consists  of  two  linear  regions,  one  extending  to 
approximately  300  %,  and  the  other,  having  a  smaller  slope,  extending  beyond  this 
thickness.  From  the  Deal-Grove  equation  the  inverse  growth  rate  is  given  by 

dt  2L  1 

dL  k  kj 

Since  k  is  proportional  to  D,  the  slope  should  be  inversely  proportional  to  D.  Fargeix 
and  cowSrkers  maintain  that  the  only  way  to  explain  the  larger  slope  for  oxide 
thicknesses  below  300  X  is  by  a  decrease  in  D  for  this  region  near  the  interface  (59,60). 
This  decreased  diffusivity  is  attributed  to  stress  in  the  oxide  which  is  a  maximum  at  the 
interface.  The  larger  diffusion  constant  for  thicker  oxides  results  from  a  relaxation  of 
the  stress  with  time.  At  any  time  the  diffusivity  in  the  oxide  is  assumed  to  have  the  form 

D(x)  =*  0^(1  -  ocexp(-x/X)) 
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where  X  is  a  damping  length  measured  from  the  outer  oxide  surface.  Using  this  form  for 
the  diffusivity  in  the  Deal-Grove  equation,  they  obtain  a  fit  to  the  experimental  data 
plotted  as  dt/dL  versus  L  by  correctly  choosing  ac  and  X  values.  Although  this  model 
claims  to  explain  the  oxidation  kinetics  in  the  initial  regime,  there  is  an  inconsistency 
in  the  logic  they  develop.  The  author's  cite  work  by  Marcus  and  Sheng  (90)  as  providing 
evidence  for  a  decreased  oxidation  rate  due  to  compressive  stress.  However,  they  claim 
this  same  compressive  stress  is  responsible  for  a  lower  diffusivity  near  the  interface 
and  thus  for  the  rapid  initial  oxidation  rate  observed  in  dry  0 ^  oxidations.  In  fact  the 
initial  assumption  of  the  model,  that  the  diffusivity  must  be  lower  near  the  interface  in 
order  to  explain  the  shape  of  the  dt/dL  plot,  maintains  that  the  Deal-Grove  model  is 
valid  in  the  initial  regime  when  it  is  known  that  this  is  precisely  the  region  where  the 
Deal-Grove  model  does  not  apply. 

It  has  not  been  clearly  resolved  whether  stress  affects  the  diffusion  process,  the 
reaction  process  or  both,  but  it  should  be  pointed  out  that  if  stress  affects  the 
oxidation  process  it  is  important  to  recognize  whether  the  reacton  is  proceeding  under 
diffusion  control  or  reaction  control  before  developing  a  model  to  explain  the  effects  of 
stress.  Clearly  experimental  values  for  stress  as  a  function  of  oxidation  temperature, 
orientation,  and  time  are  essential  for  evaluating  these  models. 


A  New  Approach 


The  most  recent  treatment  of  Si  oxidation  data  produced  a  novel  approach  to 
analyzing  characteristics  of  the  oxide  growth  kinetics  in  the  initial  regime  (49).  Data 
collected  over  the  temperature  range  800°-1000,,C  in  a  dry  0^  ambient  was  analyzed  in 
terms  of  the  rate  enhancement  in  the  initial  regime  over  that  predicted  by  the  Deal-Grove 
model.  Comparing  the  experimentally  derived  dL/dt  versus  L  plots  with  the  rates  predicted 
by  the  Deal-Grove  model,  it  was  determined  that  the  "excess"  growth  rate  could  be 
expressed  in  terms  of  two  decaying  exponentials.  In  other  words,  the  total  growth  rate  is 
the  sum  of  the  Deal-Grove  growth  rate  plus  two  additional  terms  which  decay  exponentially 
with  thickness.  One  term  decays  rapidly  with  thickness  from  the  interface  and  is  thought 
to  be  related  to  surface  impurities  possibly  introduced  by  wafer  cleaning  procedures.  An 
intermediate  phase  of  the  enhancement  is  represented  by  a  term  which  decays  more  slowly 
with  a  decay  length  that  is  roughly  independent  of  temperature.  This  behavior  was 
observed  for  all  orientations  and  temperatures  investigated.  Details  of  the 
characteristics  of  these  two  terms  as  a  function  of  temperature,  pressure,  orientation 
and  doping  level  have  been  published  (49).  It  was  pointed  out  that,  if  these  additional 
terms  arise  from  additional  fluxes,  they  must  be  in  parallel  with  the  fluxes  enumerated 
in  the  Deal-Grove  model.  Possible  sources  for  the  additional  fluxes  are  presented  in  the 
analysis.  This  approach  to  analyzing  the  initial  regime  may  provide  useful  insight  into 
those  parameters  which  enhance  the  oxidation  rate  and  should  be  incorporated  into  an 
oxidation  model. 
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A  variety  of  the  most  often  quoted  Si  oxidation  models  have  been  discussed.  Even 
many  of  these  apparently  popular  ideas  were  eliminated  from  further  consideration  because 
of  the  lack  of  agreement  with  the  presently  known  experimental  studies  of  Si  oxidation. 
Yet  all  of  these  models  agree  with  some  aspect  of  the  oxidation  phenomena  and  all  can  be 
made  to  fit  the  L, t  data.  The  danger  of  using  the  fit  of  the  film  growth  data  to  the 
model  and  too  few  facts  then  becomes  apparent.  Essentially  one  arrives  at  necessary  but 
not  sufficient  evidence.  It  should  also  be  clear  that  the  presently  known  facts  are 
insufficient  to  determine  an  unambiguous  and  complete  model  for  Si  oxidation.  It  is  then 
appropriate  to  end  this  review  with  a  brief  mention  of  some  of  the  important  missing 
parts  to  the  Si  oxidation  puzzle. 

On  the  issue  of  oxidant  transport,  it  is  necessary  to  have  a  measurement  of  the 
chemical  diffusivity,  D,  of  0^  through  SiOj  films  with  differing  film  growth  conditions 
and  various  temperatures.  At  che  present  time  only  scant  data  exists  for  bulk  fused 
silica  and  none  for  SiO_  films.  Along  with  D,  a  measurement  of  the  0.  solubility,  S,  is 
required.  No  measurements  now  exist.  From  D  and  S  the  permeability,  K,  can  be  calculated 
and  compared  with  the  values  extracted  from  L, t  data  via  various  models. 

With  regard  to  the  chemical  reaction  between  Si  and  0  at  the  Si-SiO^  interface,  the 
kinetic  order  of  the  reaction  is  unknown,  and  the  quantitative  effects  of^the  Si  surface 
orientation  and  surface  cleanliness  are  thought  to  be  crucial  parameters. 

The  role  of  intrinsic  stress,  o~ a  most  promising  new  idea,  needs  further 
elucidation  through  its  measurement  as  a  function  of  T,  P,  Si  orientation,  inert 
annealing,  SiC>2  film  thickness  and  impurity  content.  Along  with  stress  measurements  the 
measurement  of  the  viscous  relaxation  phenomenon  is  required  so  as  to  be  able  to  sort 
among  the  various  stress  models. 

The  interrelationship  of  and  to  stress,  diffusivity,  film  thickness,  growth 
and  annealing  conditions  and  surface  condition  is  rather  vague,  yet  some  relationships 
seems  to  exist  and  will  undoubtedly  contribute  to  our  understanding  of  the  Si-SiO^ 
electrical  properties.  ^ 

Finally,  the  structure  of  SiO_  films  and  the  transition  layer  between  Si  and  Si02 
will  shed  light  on  the  nature  of  tne  oxidation  process. 

Judging  from  the  recent  proliferation  of  oxidation  models,  it  appears  that  the 
current  trend  in  Si  oxidation  modelling  is  to  speculate  upon  the  oxidation  process  rather 
than  perform  the  difficult  measurements  discussed  above.  Any  experimental  progress  in 
determining  these  parameters  will  be  of  great  importance  in  resolving  the  kinetic 
mechanism  for  Si  oxidation,  a  critical  concern  of  the  semiconductor  industry. 
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